Z'^_j^ phenomenology at LHC 

Y. A. Coutinho 
Institute de Fisica-Universidade Federal do Rio de Janeiro 
Av. Athos da Silveira Ramos 149 
Rio de Janeiro - RJ, 21941-972, Brazil 

E. C. F. S. Fortes, J. C. Montero 
Institute de Fisica Teorica-Universidade Estadual Paulista 
R. Dr. Bento Teobaldo Ferraz 271 
S&o Paulo - SP, 01140-070, Brazil 

Abstract 

We study the Z' phenomenology for two extensions of the Elec- 
troweak Standard Model (SM) which have an extra U{1)b-l gauge 
factor. We show the capabilities of the LHC in distinguishing the 
signals coming from these two extensions and both of them from the 
Standard Model background. In order to compare the behavior of 
these B ~ L models we consider the reaction p + p — > /Lt+ + /z~ + X 
and compute some observables as the total cross sections, number of 
events, forward-backward asymmetry, final particle distributions like 
rapidity, transverse momentum, and dimuon invariant mass, for two 
LHC regimes: ^{C) = 1 TeV (Ifb^^ ) and 14 TeV (100 fb"^) for 
Mz' = 1000 GcV and 1500 GeV. We show that by using appropri- 
ate kinematic cuts some of the observables considered here are able 
to extract different properties of the Z' boson, and hence providing 
information about to which B — L model it belongs to. 
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1 Introduction 



The starting up of the Large Hadron Cohider (LHC) operating in an en- 
ergy range far above the electroweak scale, offers a possibihty to reveal new 
phenomena and to explore the phenomenology extracted from the expected 
huge amount of experimental data. In this way, it is expected new physics 
manifestations throughout the appearance of new degrees of freedom such 
as new charged and neutral fermions, superpartners, new gauge bosons, and 
Higgs scalar (s). To take into account these new possible degrees of freedom 
we must go beyond the SM. There are many ways to construct extensions 
of the electroweak SM. However, if we are concerned with new gauge bosons 
we must consider a larger gauge group. Whatever the extension of the SM 
is, it must present the same results of the SM at the electroweak energy 
scale in order to be consistent with present data. 

Although hadron colliders present a high background (pile-up) compared to 
electron-positron (beamstrahlung) colliders, they seem to be more suitable 
for discovering new degrees of freedom. Since high centre-of-mass energy and 
high precision can not be achieved simultaneously in the same experiment, 
both types of colliders are needed in order to have a detailed description of 
a particular phenomena. This is easily realized if we look to a recent past, 
when the synergy between hadron and lepton colliders, brought up com- 
plementary information for each other, generating benefits to high energy 
physics. As examples we have the discovery of the Z boson, the gluon, and 
the quark top. The discovery of the Z boson occurred in a "pp collider (SPS) 
but its properties (decay widths, couplings to fermions, asymmetries, mass, 
etc..) were measured with high precision in lepton colliders (SLC, LEP). 
In this way, the recent studies started at LHC can direct the search for new 
physics in future linear colliders (ILC/CLIC). 

Among the main goals of any accelerator proposal we find the physics of 
an extra neutral gauge boson Z' . This new boson is predicted by many 
extensions of the Standard Model like: 

• The models from group, also know as rank 5 models (ER5M) [HIH]; 

• The Left-Right which are based in the group SU{2)l (8> SU{2)ji (g) 

• The models with Z' in Little Higgs scenario [6]; 

• The Sequential Standard models - SSM (7j; 

• The Kaluza Klein (KK) models, predicted by extra dimensions [8]; 
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• The models with SU{3)c x SU{3)l x U{1)n gauge symmetry, known 
as 3-3-1 models [HKIOKn]; 

• Models where a strong dynamics is involved in the electroweak spon- 
taneous symmetry breaking. We can cite, for instance, Topcolor and 
BESS (Breaking Electroweak Symmetry Strongly) models |12] . 

Some models which have an extra U{1) gauge group [13^ [T^ also deserve 
a detailed study in order to predict new physics behavior. The cumulative 
study of different models, candidates to describe the electroweak interaction, 
can guide the LHC data analysis if in fact a Z' is discovered. The extra U{1) 
group gives rise to the Z' vector boson. The new constraints on these models 
generate few free parameters, if we compare them to other models which have 
origin in larger gauge groups. Following the arguments here exposed, first 
we have considered the Z' boson observables for electron-positron colliders 
|15| . and now we extend the study for hadron colliders. 
In this paper, we study the phenomenology of the Z' boson originated from 
two different extensions of the electroweak SM with an extra U{1) gauge 
factor. We predict the possible signals of new physics in the process p+p — > 
(7, Z, Z') — > fi'^+fj.~+X, at the tree level. We compare several quantities 
as total cross sections, number of events, and decay widths assuming that 
both U{1) factors correspond to the local B — L symmetry, and that the 
breaking of this symmetry occurs at an energy scale above the TeV scale. 
In order to have a more realistic scenario, we choose the parameters and 
constraints of the two different models to compare them at the same Z' 
boson mass. We use Mz' = 1500 GeV, but we also show the results for 
Mz' = 1000 GeV, for two centre-of-mass energies: ^/s = 7 and 14 TeV. 
This article is organized in the following way. In section II we introduce the 
models. In section III we present the results and in section IV we present 
our conclusions. 

2 The models 

In this paper we consider two extensions of the electroweak SM which have 
an extra U{1) local factor, resulting in the total gauge symmetry SU{2)l (8) 
U{l)i [7(1)2. More specifically, we are concerned with models which are 
based on the gauge symmetries: 

• SU{2)l O U{1)y ® U{l)z, called as Secluded model [l3]; 

• SU{2)l (S) U{1)y' <S) U{1)j,_^ called as Flipped model [H]. 
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In the Secluded model, z is a new U{1) charge, and Y is the electroweak SM 
hypercharge. In the Flipped one, Y' is a new extra U{1) charge, and the 
B — L charge assignments are those of the SM. In this last model the SM 
hypercharge is recovered after the first spontaneous symmetry breaking. 
In the Secluded model, the electric charge has no component in the U{l)z 
and in particular, depending on the value of this charge, several versions of 
this model can be studied. Our analysis will be focused on comparing both 
models only when the z charge of the Secluded model coincides with B — L. 
Besides, in this minimal version, considered here, both models have the same 
scalar sector, composed by one doublet, with hypercharge Y = +1, and a 
complex singlet with y = 0. The number of right-handed neutrinos is three, 
in order to be consistent with the anomaly cancellations. We had set them to 
be heavy enough in order to have Z' decay into heavy right-handed neutrinos 
kinematically suppressed. Besides that, the minimal version considered here 
is also absent of exotic decay channels like heavy quarks or heavy charged 
leptons. So, our results of decay widths, branchings fractions for these two 
particular Z' models are related only to standard model fermions. We will 
parameterize the neutral currents in terms of the mass eigenstate fields as 
follows: 

^^"^ = E V^7m[(5^ - 5^75)^r + {fv - fAl,)Z!^]^^. (1) 

In this paper we briefly present both models. The complete details about 
both of them were subject of study in a recent work [15]. Moreover, the 
scalar sector of this sort of models, with the respective Feynman rules, is 
considered in reference [16] . The charge operators for the Secluded and the 
Flipped models, are respectively given by: 

Q 1 

e 2 ' 



f = /3 + ^ [Y' + {B-L)]. 

When the z charge of the Secluded model is equal to i? — L it implies 
that the parameter zh (the z-charge of the SM Higgs doublet), discussed in 
Refs. \\.'6\ I14j . is equal to zero. In this case, there is no tree level mixing in 
the mass matrix between the neutral gauge bosons i?^, W!^ and By, and 
we have Z = Z\ and Z' = Z2, i.e. the mass eigenstates and the symmetry 
eigenstates coincide. This is not the case for the Flipped model. Here, we 



4 





Neutrinos 


Leptons 


u — quarks 


d — quarks 


fv 


0.8420 


0.4977 


-0.0511 


-0.3955 


fA 


-0.1732 


0.1715 


-0.1732 


0.1732 



Table 1: Z' coupling constants to the SM fermions for the Flipped model 
for Mz' = 1500 GeV. 

will assume that the mixing between the two neutral vector bosons is small 
enough so that Z k, Z\ and Z' k, Z2. 

Both models are well motivated since, in general, they can be thought as 
an intermediate symmetry, remnant of a larger (unknown) gauge symmetry, 
that was broken at a very large energy scale. Moreover, they are simpler 
compared to other class of models with an extra Z' neutral vector boson: 
simpler gauge groups and fewer free parameters. 

3 Results 

We have first chosen M^/ = 1500 GeV in order to compare the results 
coming from the two different models. For the Secluded model, this choice 
respects the bounds already established on the Z' boson mass and its cou- 
pling constants by the direct search at the Tevatron [T7] , by the electroweak 
precision tests (EWPT) at LEP II, and by the low-energy neutral current 
experiments jl8| [T9l [20] , which are well studied in the literature. On the 
other hand, the Flipped model has not received the same level of attention. 
We find that it also deserves a detailed study to establish constraints on its 
parameters and numerical couplings. 

For the Flipped model we use the following inputs: gg_]^ = 0.6132, g' = 
0.4400, u = 1987 GeV, v = 246 GeV, a = 127.9, and s^y = 0.23122 [21]. 
With these parameters we obtain Mz' = 1500 GeV, Tz' = 28.39 GeV, and 
the values of the neutral current coupling constants shown in Table 1. The 
Z couplings with SM fermions remain the same so we emphasize only the 
new contribution from the Z' boson. 

For the Secluded model we use as inputs: zh = 0, Qz = 0.2, Zg = = 1/3, 
u = 7500 GeV, and v = 246 GeV. With these parameters we obtain Mz' = 
1500 GeV, Tz' = 4.04 GeV, and the values of the neutral current coupling 
constants shown in Table 2. Due to the choice zh = in the Secluded model, 
the Z' boson has only vector couplings to fermions, and presents a leptofilic 
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Neutrinos 


Leptons 


u — quarks 


d — quarks 


fv 


0.2690 


0.2690 


-0.0897 


-0.0897 


fA 















Table 2: Z' coupling constants to the SM fermions for the Secluded model 
for Mz' = 1500 GeV. 



Mz' = 1500 GeV 


Flipped B - L 


Secluded B-L 


z' ^ J2, ^'^^ 


36% 


37.5% 


z' ^ 


18.6% 


37.5% 


z' 


42.6% 


25% 


z' w+w- 


2.8% 


0% 



Table 3: The Z' partial decay widths for fermions and charged vector boson 
for both B — L models. 

behavior around the peak region. The leptofilic character means that the Z' 
boson couples preferentially to leptons, and that the cross section for lepton 
production presents a higher peak compared to quark cross sections. This 
behavior can be seen with more evidence in + e~ — > I + 1 process, where 
/ is a lepton. 

We can realize through the partial decay widths that, as we said before, for 
the Secluded model, Z' couples preferentially to leptons, as shown in Table 
3. Due to this character, and depending on the choice of the Z' parameters, 
the Secluded model can give an explanation for the positron excess in cosmic 
rays, as presented by the PAMELA experiment |22j . 

We present in Figure 1 the total Z' decay width against Mz' ■ Both of them 
are linear functions of Mz', as it should be, and hence they double their 
values in the Z' boson mass range showed in the figure. However, according 
to the parameter choice, gz = 0.2, the decay width of the Secluded model is 
smaller than that of the Flipped one. In the latter model the Z' boson has 
both axial and vector couplings, and in the Secluded model it has only the 
vector ones due to the choice zh = 0. 

The Z' production can be achieved by the Drell-Yan mechanism which im- 
plies in leptonic decay products accompanied by background, which con- 
taminates the observables under consideration. This background can be 
minimized applying some specific cuts. To perform our calculations we 
use the package Comphep [23] with Cteq6Ll parton distribution functions. 
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Figure 1: Evolution of total decay width with Z' boson mass. 

For the first simulations we have used the following cuts on the final lep- 
tons, which we call as the first cut set: pt^ > 20 GeV, |r/^| < 2.5 and 
-0.99 < cos 6*5^ < 0.99. 

The first observable obtained, when the first cut set is used, is the total cross 
section for p + p — )■ (7, Z, Z') — [i^ ^ pT ^ X. Considering Mz' = 1500 
GeV we obtain for the Flipped, Secluded and SM, 1.33 x 10^ pb, 1.01 x 10^ 
pb, and 0.6 x 10^ pb, respectively for ^/s = 14 TeV, which implies consider- 
ing the annual luminosity £ = 100 fb~^ in a huge number of events: ~ 10^. 
We present in Table 4 the number of events considering Mz' = 1500 and 
1000 GeV for both models and also y/s = 14 and 7 TeV. The number of 
events will decrease when we apply more restrictive cuts. In order to assign 
the correct quark direction, we had selected dimuon large rapidity events, 
lUfifil > 0.8, as suggested in [24l [25], the other cuts are ptfi > 20 GeV, 
|r?^| < 2.5, 1450 GeV < M^^ < 1550 GeV, and -0.99 < cos6lg^ < 0.99, and 
we call them as the second cut set. 

Besides the total cross section and decay width for leptonic final state, we 
present: invariant mass, rapidity, forward-backward asymmetry, transverse 
momentum and angular distribution. In Figures 2a and 2b we present the 
dilepton invariant mass distribution. The Flipped model presents a distinct 
behavior with a higher peak far above the Standard model and has very 
good changes to be separated from the background. On the other, the 
Secluded model distribution in the region off-peak in Figure 2 (down panel) 
has a behavior similar to that of the SM, while the on-peak region is not 
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Model 


Mz' = 1500 GeV 


Mz' = 1000 GeV 


Mz' = 1000 GeV 


Flipped B - L 


1.33 X 10** 


8.48 X 10^ 


245 


Secluded B -L 


1.01 X 10" 


7.18 X 10'' 


219 



Table 4: Number of events considering Mz' = 1500 and 1000 GeV, and 
the first cut set. The second and third columns are for y/s = 14 and C = 
100 fb~^, and the forth column is for ^/s = 7 TeV and C = Ifb^^. In the 
data related to Mz' = 1000 GeV we have also set M^^ > 500 GeV. 

so powerful to distinguish it from the SM. This difference has origin mainly 
in the choice of zh = for the Secluded model, which leads to pure vector 
couplings of the Z' boson to fermions. A different choice for zh, say zh ^ 0, 
will restore the axial couplings but, in this case, the z-charge will not be 
equal to B — L. The case of z// 7^ is out of the scope of the present work. 
The rapidity distributions can also be used to disentangle the B — L models. 
This observable is powerful even when we consider both light and heavy 
quarks or only the light ones, this result is showed in Figures 3a, 3b and 3c. 
The signal has been improved for both models but once again the Flipped 
model produces a signal above the ones of the other two models. One in- 
trinsic characteristic of the Flipped model is the increasing signal until it 
reaches the region where the rapidity of each muon is \y^\ = 1.25. Both, the 
Secluded model and the SM do not have this behavior. When we apply more 
severe cuts on rapidity, as |?/^| > 0.8, the resulting data reveal that the three 
models have less possibilities to be disentangled. So, for this observable, the 
first cut set is more efficient than the more severe ones. 
Another important observable, very sensitive to the new physics contribu- 
tions, is the forward-backward asymmetry {A^p^). In e^e~ colliders, A^p^ 
is measured with high precision, due to the well known initial beam direc- 
tion. However, the situation is quite different for hadron colliders where the 
original quark direction is completely unknown. To solve this problem, the 
quark direction can be approximated by the boost direction, which connects 
the dimuon reference system with the original beam direction. In order to 
achieve this goal we have applied the \y^^\ > 0.8 cut for the muon pair ra- 
pidity, following the references |24tl25|. In contrast with the Flipped model, 
the Secluded one has no Z' axial couplings to fermions and, as a consequence 
of it, the forward-backward asymmetry, even being a powerful tool to get 
information on the Z' boson parameters, does not receive contributions from 
the axial couplings on-peak. In the Figure 4 we show this observable versus 
M^^, and we see that the two B — L models can be distinguished. 
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If we consider all the studied observables, after applying the second cut set, 
we will realize that one particular numerical cut can put in evidence one 
observable and others not. The cut > 0.8 is not so helpful for rapidity 
distributions, but it is necessary for forward-backward asymmetries in order 
to guarantee that the quark direction has more probability to be the correct 
one. It means that the probability of obtaining more significant measure- 
ments of forward-backward asymmetries increases for higher rapidities. 
We present in the Figure 5 our results for the muon transverse momentum 
distribution. One of the characteristics of this distribution is that the peak 
is located right at the position px = Mz' /^^ for both B — L models. We 
have adopted more restricted cuts in the final dimuon mass, M^^ > 500 
GeV [26] in order to emphasize this behavior. The flipped has once again, 
much better chances to be separated from the background, in opposition, 
the Secluded one can be mistaken as Standard model background. 
The muon angular distribution is less sensitive to Z' contributions, but can 
be used as a previous result before the calculation of A^p^. From Figure 
6 we can realize that the signals related to the SM are below the signals 
of the two B — L models. We can see a small asymmetry when the curves 
cross each other in the region —0.50 < cos0 < —0.25, where 9 is the angle 
between the beam and the muon directions. 

We now consider a Z' boson which mass is 1000 GeV and y/s = 7 and 
14 TeV in order to extract more conclusions about the same observables 
presented above. The related graphs are presented in the Figures 7, 8, 9 
and 10. The increase of the centre of mass energy from 7 to 14 TeV has as 
direct consequences the increase of the cross sections and, hence, the number 
of events. The invariant mass distribution shows also a small peak for the 
Secluded model if compared to the Flipped one, as we can see in Figures 7a 
and 7b. The rapidity distributions presented in Figures 8a and 8b show that 
for Mz' = 1000 GeV both models can not be clearly distinguished from each 
other in this case. If we consider = 14 TeV the rapidity distributions 
graphs are separated, mainly in the region < y^^ < 1, but this separation 
is not so noticeable as the one for Mz' = 1500 GeV, already presented in 
Figures 3a and 3b. The transverse momentum distributions, presented in 
the Figures 9a and 9b, do not show a clear peak for Mz> /2 for the Secluded 
model, but the Flipped one is put in evidence once again with higher values 
for this distribution on-peak. The muon angular distributions present a 
better behavior for disentangling both models, mainly if we consider ^/s = 14 
TeV, as can be seen from Figures 10a and 10b. 



9 



4 Conclusions 



In order to explore the properties of a new neutral vector boson belonging to 
models with B — L local symmetry, we study the Drell-Yan channel with a 
muon pair production, for two LHC proposed scenarios: ^/s = 7 TeV (£ = 1 
fb~^) and y/s = 14 TeV (£ = 100 fb~^). Our study is focused on two values 
of the Z' boson mass, 1000 GeV and 1500 GeV. 

First we can say that if Mz' is smaller than the centre of mass energy of 
the colliding protons, it can be discovered at the LHC. For instance, for 
Mz' = 1500 GeV and ^/s = I'i TeV we obtain ~ 10® events for both B - L 
models. The signals of Z' bosons come from their direct production. 
In order to perform Z' studies in pp collisions, many observables can be used. 
We can cite the total cross-section and also several asymmetries very well 
studied by e^e~ colliders. The measurement of some of these observables 
is very difficult due to the signal, which needs to be detected above the 
background of the hadronic experiment. Thus, it is necessary to apply cuts 
taking into account the detector acceptance for reducing the background 
contribution. 

The invariant mass distribution of the muon pair final state is centered 
around the Z' mass. This behavior allows us to separate the signal from the 
background and to measure the Z' boson mass. Several cuts can be applied 
to guarantee more precision on measurements. The signal from the Flipped 
model is above the one from the Secluded model, and the corresponding 
difference can be enhanced if more severe cuts are considered. On the other 
hand, the SM background remains below the signals of both B — L models, 
as shown in Figure 2. 

The muon transverse momentum distribution centered around the half of 
Z' boson mass is a helpful observable for distinguishing the models. In this 
case we have observed a dominance of the Flipped model over the Secluded 
one, which has a behavior similar to the SM. See Figure 5. 
The forward-backward charge asymmetry is another powerful tool to provide 
some information about the Z' couplings to quarks in pp collisions, when 
leptons are in the final state, even being measured indirectly due to the 
unknown of original quark direction. To guarantee the correct direction of 
the initial parton we applied a cut on final muon pair > 0.8. The 

resulting ApB asymmetries, for both B-L models, present clear different 
behaviors, mainly for the muon pair invariant mass in the range 1300 — 1500 
GeV, as shown in Figure 4. 

In addition to ApB, the rapidity distributions are very sensitive to the Z' 
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couplings to quarks and in their computing we can consider light and heavy 
quarks. They are very useful to disentangle the different models. Different 
cuts on it should be analyzed to exploit different responses from the same 
experiment. The combination of the forward-backward asymmetry and ra- 
pidity distribution generates another interesting observable, Ap^yi- If the 
LHC discovers a Z' boson, the ApByi can be used as a refinement to select 
to which model the Z' boson belongs to. This combined observable in not 
under consideration in this work. 

Besides the study for LHC operating in the high energy regime, we have 
performed an analysis at ^/s = 7 TeV for AIz' = 1000 GeV. The obtained 
results present a clear signature for the Z' boson production. The dimuon 
invariant mass distribution shows a sharp peak around Z' boson mass, above 
the SM background, as shown in Figure 7. We note that the transverse 
momentum distribution once again reveals a peak on pT = Mzi/2. For 
these two distributions, the ones of the Flipped model are far above those 
of the Secluded model. 

In general, the presented models can be well exploited by using the cited 
observables above at the LHC first energy stage. Some of them are more 
suitable for distinguishing the models. The choice of the applied cuts is also 
important. For instance, the same cuts that we had used to the forward- 
backward asymmetries are not appropriate to the rapidity distributions. 
Therefore, a combined analysis and different cuts should be applied to dis- 
entangle the models and guarantee more realist theoretical prediction on the 
Z' couplings to fermions. The LHC running at y/s = 7 TeV already is a 
powerful tool for unraveling physics beyond the SM. 
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Figure 2: The dilepton invariant mass distribution for p + p — ?• + +X 
process for the Flipped, Secluded and Standard models considering Mz' = 
1500 GeV at ^/s = 14 TeV. First cut set (up) and second cut set (down). 



14 



4501 ' 1 ' ^ ' 1 ' ^ ' 1 ' r- 




M^. = 1500 GeV 


B-L Flipped 

B-L Secluded 

SM 




'v. 












^' s = 14TeV 





50 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' '— ^ ' 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

y valence 



p + p^H+ + n +x 



10^ 



10' 



M2. = 1500 GeV 


B-L Flipped 




B-L Secluded 




SM 




Vs = 14TeV 







0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

y 



Figure 3: The rapidity distribution for the process p + p— )•//"'' + + X 
considering Mz' = 1500 GeV for the Flipped, Secluded and Standard models 
at ^/s = 14 TeV. The valence and sea quark contributions (up) and only 
the valence quarks (middle), using the first cut set. The same rapidity 
distribution using the second cut set (down). 
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Figure 4: The forward-backward asymmetry for the process p -\- p ^ + 
fi~ + X as a function of the dilepton invariant mass for both B — L models, 
considering Mz' = 1500 GeV, ^/s = 14 TeV, and applying the second cut 
set. 
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Figure 5: The muon transverse momentum distribution for the process p + 
p — >• ^'^ +n,~ +X for the FHpped, Secluded and Standard models, considering 
Mz' = 1500 GeV, ^/s = 14 TeV, and applying the second cut set. 
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Figure 6: The muon angular distribution for the process p+p — )• /i + +X 
for the Fhpped, Secluded and Standard models, considering Adz' = 1000 
GeV and = 14 TeV. 
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Figure 7: The dilepton invariant mass distribution for p+p — ?■ + fi~ + X 
process in Flipped and Secluded models considering Mz' = 1000 GeV at 
^/s = 7 TeV (up) and ^/s = 14 TeV (down) with the first cuts set. 
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Figure 8: The rapidity distribution for the process p + p — ?■ /i"*" + fi~ + X 
for both B — L models, considering Mz' = 1000 GeV, ^/s = 7 TeV (up) and 
^/s = 14 TeV (down) with the first cut set. 
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ure 9: The muon transverse momentum distribution for the process 
p — )■ fi'^ + n~ + X ioi both B — L models, considering Mz' = 1000, 
= 7 TeV (up) and ^/s = 14 TeV (down), applying the first cut set. 
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Figure 10: The muon angular distribution for the process p + p ^ fi 

+ X for both B-L models, considering M^/ = 1000, ^/s = 7 TeV 
and = 14 TeV (down), applying the first cut set. 
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